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^ ■ Abstract 

We consider top-anti-top production near threshold in e + e~ collisions, resumming 
Coulomb-enhanced corrections at next-to-next-to-leading order (NNLO). We also 
sum potentially large logarithms of the small top quark velocity at the next-to- 
leading logarithmic level using the renormalization group. The NNLO correction 
to the cross section is large, and it leads to a significant modification of the peak 
position and normalization. We demonstrate that an accurate top quark mass 
determination is feasible if one abandons the conventional pole mass scheme and if 
one uses a subtracted potential and the corresponding mass definition. Significant 
uncertainties in the normalization of the ti cross section, however, remain. 



1. Introduction. The top quark mass is now known to be around 175 GeV with 
an accuracy of 5 GeV from the direct measurement at the Fermilab Tevatron Collider 
]l| . Accurate mass determinations (with errors below 1-2 GeV) are difficult at hadron 
colliders. Despite the fact that orders of magnitudes more top quarks will be produced 
at the CERN Large Hadron Collider, a precision measurement is reserved for a future 
lepton collider. In this case the method of choice relies on scanning the top quark pair 
production threshold. From an experimental point of view, an error in the 100 MeV 
range is conceivable || ; the limiting factor is the accuracy to which the cross section can 
be predicted theoretically as a function of the top quark mass. 

The literature on top quark physics near threshold in e + e~ collisions is substantial 
||. ||. Perturbative calculations in the threshold region require that either the toponium 
Rydberg energy scale m t a 2 s S> Aq CD or that the top quark decay width Y t ^> Aq CD . Both 
conditions are satisfied and since T t ~ m t a 2 s ~ 1.5 GeV, narrow toponium resonances do 
not exist ||. In the kinematic region of interest, the top quarks are slow, with typical 
velocities v ~ 1/10. As a consequence methods familiar from non-relativistic bound 
state calculations in QED can be used to compute the ti cross section. In particular, the 
dominant interaction between the t and i can be described by the colour-singlet Coulomb 
potential, which has to be treated non-perturbatively [||]. Corrections are computed in 
the background of this strong Coulomb interaction. 

Recently, the 2-loop QCD correction to the Coulomb potential (correcting an 
earlier result 0) and the 2-loop relativistic correction to the ti vector coupling to the 
initial virtual photon or Z boson || |9| have been computed. With these two inputs at 
hand, we can compute the next-to-next-to-leading order (NNLO) QCD correction to the 
top quark pair production cross section. (The precise meaning of 'NNLO' in the present 
context is given below.) In the following, we first compute the NNLO QCD correction 
in the conventional pole mass scheme, i.e. the threshold cross section is expressed in 
terms of the top quark pole mass. Comparable calculations, with some technical and 



implementational differences, have already been completed by several groups |TI| |IT|, |l2j] . 
We find that the NNLO correction is substantial and leads to an uncomfortably large shift 
in the top quark mass, raising questions as to the stability of the theoretical prediction. In 



, one of us suggested that such shifts could occur as a consequence of on-shell mass 
renormalization, which is particularly (and in the case at hand, artificially) sensitive 
to non-perturbative, long-distance effects. This sensitivity to non-perturbative effects 
can be removed by using a different mass renormalization scheme (called the potential 
subtraction (PS) scheme in [|IB|]) with the additional benefit that the new mass definition 
can also be related more reliably to short- distance masses (such as the MS mass), which 
are ultimately of more interest for high-energy processes and Yukawa coupling unification 
relations (should they exist). 

The main result of this letter is to demonstrate that this procedure works. We show 
that the mass shifts become small in the PS scheme, and that the PS mass can also be 
accurately related to the MS mass. However, the normalization uncertainty of the cross 
section remains large, and is larger than could have been anticipated from the earlier 
next-to-leading order (NLO) calculations. The second improvement which we suggest in 
this letter is to use renormalization group equations to sum large logarithms, lnt> and 
In v 2 to all orders in perturbation theory. The consideration of logarithms is important to 
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identify correctly the scales of the various subprocesses, but not restricted to the scale of 
the QCD coupling only. In the present work, we complete this program (almost trivially) 
to the next-to-leading logarithmic (NLL) order. The step to NNLL is substantially more 
complicated and we will discuss it, together with the details of the present calculation, 
in a future publication . 



Before proceeding, let us emphasize that 'NLO', 'NNLO' etc. do not refer to a 
conventional loop expansion in a s , because the Coulomb interaction cannot be treated 
perturbatively in the threshold region. To be precise, a LO, NLO, ... approximation to 
the normalized cross section takes into account all terms of order 



R = crtt/a^ =«£( — ) ' {l (LO); a s , v (NLO); a 2 , a s v,v 2 (NNLO); . . .} , (1) 
k ^ v ' 

where logarithms of v are suppressed. The renormalization group improved treatment 
extends this to the summation of logarithms such that a LL, NLL, ... approximation 
accounts for all terms of order 

R = vY,(— ) fc E(^ ln ^)'-{ 1 (LL); a s ,v(NLL); a 2 , a s v, v 2 (NNLL); . . .} . (2) 

k ^ v ' i 

The result discussed here is complete at 'NNLO' and 'NLL'. Furthermore, near the 
would-be toponium poles another resummation is necessary, which we discuss below. 

We also emphasize that the NNLO QCD correction is defined as above in the limit 
T t <C m t a 2 . When T t ~ m t a 2 , as we expect, further corrections arise from so-called non- 
factorizable contributions ]TS|, which have not been calculated so far with the Coulomb 
interaction treated non-perturbatively. Electroweak corrections also enter and the entire 
concept of a tt cross section has to be revised, since single resonant contributions are 
expected to contribute at NNLO (in the above power counting scheme with Y t ~ m t a 2 ) 
to the WWbb final state. These are interesting issues to be studied, but we expect them 
to bear little on the issue of an accurate top mass determination which we address in this 
letter. Hence, we will keep the top quark width only in the form of an imaginary mass 
term iT^i/j for the non-relativistic top quarks; this amounts to evaluating the tt Green 
function at energy E + iT t as familiar from LO and NLO calculations |3|, |J. As a final 
(trivial) simplification, we restrict our attention to ti pairs produced through the vector 
coupling to a virtual photon. The vector coupling to a Z boson can be accommodated 
by a trivial replacement of the electric charge. The axial vector coupling is suppressed 
by a factor of v 2 near threshold; no QCD corrections to it are needed at NNLO. 

2. Derivation of the cross section. With the treatment of the top quark width as 
specified above, and neglecting the axial- vector coupling, the tt production cross section 
is obtained from the correlation function 

IW) = (q^v - q%„) Tl(q 2 ) = i J Ae^ 1 <0|T(^(x)j,,(0))|0>, (3) 

where j^{x) = \t^H](x) is the top quark vector current and s = q 2 the centre-of-mass 
energy squared. Defining the usual i?-ratio R = o t il &o (°"o — ^ na l m / (3s) , where a em is 
the electromagnetic coupling at the scale 2m t ), the relation is 



Attp 2 

r = Z^i m n M ( s + ie), (4) 
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where e< = 2/3 is the top quark electric charge in units of the positron charge. Only the 
spatial components of the currents contribute up to NNLO. In the following, m t refers 
to the top quark pole mass. 

According to (|I|), at NNLO, we have to extract, to all orders in a s , the first three terms 
of the expansion of any Feynman diagram in v — ((\/s — 2mt)/m t ) 1 ^ 2 . The expansion in 
v is constructed [|16|] by dividing the loop integral into terms related to hard (Iq ~ m t , 
I ~ m t - referring to a frame where q — 0), soft (Z ~ m t v, I ~ m t v), potential (Iq ~ m t v 2 , 
I ~ m t v) and ultrasoft (Z ~ m t v 2 , I ~ m t v 2 ) momentum. This split-up requires a 
regularization to deal with divergent integrals that appear in intermediate expressions 
and we use dimensional regularization with MS subtractions. This procedure has already 
been used to obtain the cross section at threshold up to order a 2 . By expanding the 
all-order result presented below up to order a 2 , we reproduced the result of f9j, which 
has been used as a common input to the previous |T(], [11], [12] NNLO ti cross section 
calculations. 

We begin with integrating out the hard modes, commonly termed 'relativistic cor- 
rections'. The effective 7*t£ coupling seen by the non-relativistic quarks (described by 
two-spinor fields ip for t and % f° r t) after integrating out the hard modes is given by 

iff = Cl - A ^o\iD) 2 x + • • • , (5) 

where the ellipsis refers to terms not needed for H(q 2 ) and at NNLO. At NNLO, we can 
use c 2 = 1, while C\ is needed at order a 2 . The expression for c\ is 



Ci(/i) = 1 + 



a s (ii) 



2 



+ C i7f^ + --- ( 6 ) 



4tt 1 (4tt) 



where 5i, related to the 2- loop anomalous dimension of the non-relativistic current xp'^x, 
is given by 5i = — 5607r 2 /(276 ) (&o = 11 — 2rif/3, rif = 5) and C\ at the scale /i/j, at 
which QCD is matched onto a non-relativistic effective theory, is given by 



c x {Hh) = 1 ^ — + 



' 2 , 35tt 2 
-b + 




3 27 

89 511tt 2 14tt 2 , n 125C(3) 11 

In 2 H n f 

54 324 9 9 27 ; 



VP 2 ; + . . . . (7) 
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The first order result is well known ilTfl; the second order contribution is from 



Eq. (|g) follows from solving the renormalization group equation for c\ with the 2-loop 
anomalous dimension. When evaluated at a scale fi of order mv, this expression sums 
next-to- leading logarithms of the form a s (a s lnv) 1 to all orders. This is in fact the only 
source of next-to-leading logarithms in the problem (there are no leading logarithms), 
and @ is sufficient to obtain a NLL approximation. There is an ambiguity in the scale 
of a 2 in (||). In our numerics, we actually choose the expression for the NNLL- improved 
coefficient function, setting the (unknown) 3-loop anomalous dimension of the current 
to zero. 

The correlation function ([3]) is now expressed in terms of the effective current (|J). 
(This leaves out a hard correction from the region x ~ l/m t in the integral @. However, 
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in this region the top quarks are far off-shell and no contribution to the imaginary part 
of II is obtained.) The correlation functions of the non-relativistic current are then 
computed with the non-relativistic effective Lagrangian. It is a straightforward matter 
of counting powers of velocity (using the momentum scaling rules given above) to show 
that the following terms in the effective Lagrangian are sufficient at NNLO: 

£nrqcd = (iD° + ^- + it\ V + ^3 V^V -^tfv B^j 

+ ^ ^ \d\ E'] if; + -f^- ^a 11 \D\ E j ] if; + antiquark terms + £ light . (8) 



5m; 



■smi 



Because we use dimensional regularization, some care is needed to define the algebra 
of Pauli matrices in 3 — 2e dimensions as well as anti-symmetric products. We use 
= [a\ <ji}/(2i) (equal to e l ^ k a k in three dimensions) and cr ■ B must be interpreted as 
in terms of the gluon field strength tensor. The last term in (|), £n g ht, denotes 
the QCD Lagrangian of the massless fields, i.e. the QCD Lagrangian with the top quark 
part omitted. The coefficient functions di, o^, d% can be set to 1 at NNLO. Their leading 
logarithmic renormalization would be required for a NNLL approximation. In that case, 
further operators, notably four-fermion operators (of heavy-heavy and heavy-light type) 
would have to be added to the effective Lagrangian. We discuss this extension in [|j~4}| . 
The unconventional term involving the top quark width accounts for the fast top quark 
decay as discussed in the introduction. (We should emphasize again that the Lagrangian 
is not complete to NNLO as far as the treatment of the width is concerned. The term 
we added is the leading order term, but further terms exist which are suppressed by two 
powers of velocity.) 

The loop integrals constructed with the non-relativistic Lagrangian still contain soft, 
potential and ultrasoft modes. Near threshold, where energies are of order m t v 2 , only 
potential top quarks and ultrasoft gluons (light quarks) can appear as external lines of 
a physical scattering amplitude. Hence, we integrate out soft gluons and quarks and 
potential gluons (light quarks) and construct the effective Lagrangian for the potential 
top quarks and ultrasoft gluons (light quarks). Because the modes that are integrated 
out have large energy but not large momentum compared to the modes we keep, the 
resulting Lagrangian contains instantaneous, but spatially non-local interactions. In the 
simplest case, these reduce to what is commonly called the 'heavy quark potential'. We 
refer to this theory as potential non-relativistic QCD (PNRQCD), adapting the term 
PNRQED introduced in @ to the QCD case. 

The derivation of the potentials in the framework of the threshold expansion |16|] will 
be presented elsewhere [14|. The following result has been obtained by matching the 



on-shell tt scattering amplitude in NRQCD onto its PNRQCD counterpart. We verified 
that the potential is gauge-independent for a general covariant gauge and the Coulomb 
gauge. (To obtain a gauge invariant result one has to combine the contribution from 
the soft modes with the one from potential gluons.) The resulting momentum space 
potential, after carrying out a colour and spin projection on the components relevant to 
the calculation of IT, is 

V(p, q ) = -^^ + 5V(p,q), (9) 
Q 
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5V(p,q) 



(2 2 
a 2 - {2a x b Q + h) ln^ + b 2 In 2 ^ 



a: 



+ 



7ra s |g| / q 2 e lE 



4m+ 



^ J (4tt) 2 

) vrWa - 2e) ^-2-(l-2e) + C x (l-e) 



~ 2 r d 2 -7d+10 , 2 



+ — — h — 

mi 



m\ 



A{d-1) 



4 



(10) 



where d = 4 — 2e, CV = 4/3, C*a = 3, d\ = d 2 = 1 in the present NNLO/NLL approxi- 
mation, and bi = 102 — 38n//3 the two- loop coefficient of the QCD /3-function. (Always 
a s = a s (fi).) The loop corrections to the Coulomb potential are a\ = (31Ca/9 — 10n//9) 
and § 



a 2 = C 2 



4343 22C 3 2 vr 4 

I + 4tt 2 

162 3 4 



- C A n f 



81 



28C(3) 



CpUf 



T ~ 8C(3) 



100n 2 f 



(11) 



The potential ( [T0| ) differs from the potential used in []TD|, [U], |12"fl . Firstly, we need the 
potential in d space-time dimensions, because the terms in the last two lines generate 
ultraviolet divergent integrals, which we regularize dimensionally. The divergences cause 
a factorization scale dependence which cancels with the factorization scale dependence 
in the coefficient function (|B|) of the non-relativistic current. (The Coulomb potential 
does not generate divergences and we can use a\ )2 in four dimensions.) Secondly, our 
potential contains a C F a 2 s /{m t \q\) term, which is absent in |TI| [TTJ, [T2|] . Nevertheless, 
both potentials (in four dimensions) are in fact equivalent |[14|| . 

Having integrated out soft modes and potential gluons, the correlation functions of 
non-relativistic currents are computed with the Lagrangian 



pnrqcd = V> f (id + + ir t j ip + X ] (id 



2mi 



+ ir t )x 



I d d '\ 
1 



C F a s 



x ] x 



(0) 



+ 



5m, 



X^d\ + / d d ~ l r Utyl ( r ) SV ( r ) 



X X X 



(0) 



(12) 



The terms in the last line are treated as perturbations. However, velocity power counting 
reproduces the well known fact that the leading order Coulomb potential in the second 
line is not suppressed compared to the free non-relativistic Lagrangian in the first line. 
Consequently perturbation theory in PNRQCD implies that instead of freely propagat- 
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ing, a tt pair propagates according to the Coulomb Green function, which satisfies 



E] G c (p,p'; E) + 



d d ~ l k 



(2 



7T 



, d — 1 



k 1 



G c (p-k,p';E) = (2n) d - 1 6^ d - 1 \p-p') 

(13) 

with E = E + iT t and E = ^/i — 2m t . Because we use dimensional regularization, 
all quantities are defined a priori in momentum space; the above equation defines the 
Coulomb Green function in d dimensions. The PNRQCD Lagrangian (|T2|) does not con- 



tain any gluon fields any more. This is so, because at NNLO the top quarks interact only 
through potentials. Counting powers of velocity for the leading ultrasoft interactions, we 
find that they are of NNNLO and higher order, i.e. beyond the accuracy of the present 
calculation. 

To complete the calculation, we compute with PNRQCD the correlation functions of 
the non-relativistic currents. For the power-suppressed term in (||) the LO Lagrangian 
suffices. For the current \ 0"V we need the kinetic energy correction to first order, the 
potentials suppressed by one power of a s or v relative to the Coulomb potential to second 
order, and the potential suppressed by two powers of a s or v to first order. The explicit 
result for the cross section is lengthy and will be given in |14]. We have checked this 
result by expanding it to order a 2 , confirming the result of in this limit. This gives us 
confidence that the factorization in dimensional regularization has been done correctly. 
The Coulomb Green function contains bound state poles at E^° = —m t {CFCi s ) 2 /{in 2 ) 
for positive integer n. The location and residues of the bound state poles are modified 
by QCD corrections. We computed the bound state energies to order and residues to 
order oP s and find agreement with the results of and [21j], respectively. 

The calculation described so far sums correctly all terms at NNLO and NLL, as 
defined in (p]) and (@). However, the result contains terms of the schematic form 



a„E. 



LO 



E LO 

n 



(E + iT t 



(14) 



which become large in the vicinity of E = E^°, if T t is not much larger than E^°. 
For top quarks T t ~ E\° and it is necessary to sum singular terms of the form (14) 
to all orders. This is easily done by adding the expression with the exact bound state 
energy denominator at NNLO and by subtracting the same expression but expanded and 
truncated at NNLO. That is, we add: 



(1 + f ias + f 2 a 2 ) 

£LO (x + eias + e2a 2) _ e \E^°-E 



- Ef E^ _ E 



F^{E^fe\ F n LO E n LO (e 2 + e 1 / 1 ) + F^f 2 



(2#° - Ef 



(£LO _ Ef 



El: - E 



(15) 



This procedure (also discussed in |TT|, P2| ) is essentially equivalent (near the bound 
state poles) to solving the Schrodinger equation exactly with the potential SV(r), rather 
than treating it perturbatively as we have done so far. In practice, we implement this 
resummation for the first two bound state poles. The correction for the remaining ones 
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is tiny, because the residues decrease as 1/n 3 . It is worth noting that even after this 
resummation, one would not be able to compute the ti cross section in the threshold 
region, if the width of the top quark were smaller than about two times the error that 
remains in the location of the IS bound state pole at NNLO. In the conventional pole 
mass scheme this requires T t to be larger than roughly 1 GeV, a constraint which is 
satisfied but not by a large margin. 

3. Top quark mass definitions, the PS scheme. The top quark cross section 
at LO, NLO and NNLO (including the summation of logarithms at NLL) is shown 
in Fig. [I|a. (To be precise, the NLO curves include the second iterations of the NLO 
potentials.) The NNLO correction is seen to modify the line shape at the level of 20%. 
It also shifts the position of the peak by approximately 600 MeV. This conclusion is in 
qualitative agreement with the results of |L0], |TTJ, [L2]]. We should note, however, that 
our result is implemented in a different way: for instance, we do not keep the short- 
distance coefficient Ci(/x) 2 as an overall factor, but multiply it out, keeping all terms to 
NNLO. Furthermore, we prefer to choose a different renormalization scale (equivalent to 
the 'soft scale' in |10|, [TJJ |12[), typically of order 30 GeV, compared to the central value 
75 GeV chosen in previous works. This choice is motivated by the fact that the typical 
momentum transfer in the instantaneous interactions is of this order, or, if anything, 
smaller. 

The significant shift in the location of the peak impacts directly on the top quark 
mass measurement. There is no unique choice of the concept of the top quark mass. The 
top quark pole mass definition has been universally assumed in previous cross section 
calculations near the threshold; this is indeed an intuitively plausible choice as top quarks 
do not hadronize. However, the top quark pole mass definition is known to be more 
sensitive to non-perturbative effects P3| than other mass definitions and the finite width 
of the top makes no difference in this respect. (The difference is, that the top quark pole 

In 



mass is irrelevant, because a top quark is always off-shell by an amount y/m t T t .). In [TT3J 
we argued that the shift of the peak position is related to large perturbative corrections 
which appear only, when the cross section is expressed in terms of the pole mass, and 
which have their origin in the sensitivity to distances larger than the toponium Bohr 
radius. The point is that the Coulomb potential in coordinate space and the pole mass 
receive the same large corrections ||T3"| , |24jj . We therefore perform a subtraction on the 
potential such that the subtracted terms are absorbed into a mass redefinition and at 
the same time cancel the large corrections to the pole mass. This leads to the 'potential 
subtraction (PS) scheme' and the corresponding mass definition. For further details of 
the argument we refer to [13]. The PS mass at the subtraction scale /i/ is defined by 



m ti ps(At/) =m t - Sm t (iif), 



(16) 



where 

8m t {yL f ) 



+ 



1 
2 



d 3 q 

(27P 



[V(q)] 



Coulomb 



A(/V) 



o:. s 



a 2 - {2ai& + M hi 



4-2) +6g (ln 2 4-41n4 + 
fi 2 I V fi 2 fi 2 



(17) 
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Figure 1: (a) [upper panel]: The normalized it cross section (virtual photon contribution 
only) in LO (short-dashed), NLO (short-long-dashed) and NNLO (solid) as function of E = 
\fs — 2mt (pole mass scheme). Parameters: m t = fj-h = 175 GeV, T t = 1.40 GeV, a s (mz) = 
0.118. The three curves for each case refer to fi = {15(upper); 30(central); 60(lower)} GeV. 
(b) [lower panel]: As in (a), but in the PS mass scheme with [if = 20 GeV. Hence E = 
yfs — 2m ti ps(20 GeV). Other parameters as above with m t — > m 4i ps(20 GeV). 
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and A(/i/) = Cpoi s iJLf/'K. In the following, we re-express the tt cross section in terms of 
the PS mass. We suppose that \Xf scales as m t a s and count A(/ij) as order m t v 2 . Then all 
terms are re-expanded and terms beyond NNLO are dropped (modulo the resummation 
near the bound state poles discussed above). Note that A(/i/) is not expanded when it 
occurs in the combination s/s — 2(m t) ps(/x/) + A(/i/)). It should also not dominate this 
combination and this is why we do not use the MS mass directly, which would lead to 
A ~ m t a s . 

The peak in the it cross section profile is, roughly speaking, the remnant of the first 
bound state pole (the 'IS pole'). To understand the effect of the mass redefinition on the 
cross section qualitatively, it is useful to compute the correction to the IS pole in terms 
of both mass definitions. The dominant correction in the pole mass scheme is related to 
the running coupling in the LO Coulomb potential. So let us take 

rrV / x -4:iiC F a s f b a s q 2 \ k 

to compute the energy level shift 

SE k = I .f* ,f* n»(P + <z/2) SV h (q) *os(p - 9/2) 




Q 2 



{m t C F a s ) 2 



5V k {q). (19) 



The mass of the IS bound state, including only the leading order Coulomb interaction 
and the above perturbation, is given by 



M 1S = 2m t + E\° + 5E k = 2m ttPS (fi f ) + El°> PS (n f ) + [5E k + 25m t (fi f )] . (20) 

Comparing ( |19"D with (]17D , we see that the integrands in 5E k + 28m t {iXf) cancel each 
other for \q\ < ra t Cpa s and \q\ < fif. At the same time the integral (|T9| ) is dominated 
by the contribution from small q quickly as k increases. The result of this exercise for 
H = 30 GeV (in which case m t CFa s (ii) is also about 30 GeV) and \ij = 20 GeV is shown 
in the following table: 



k 


1 


2 


3 


4 


5 


6 


7 


8 


5E k /MeV 


-489 


-214 


-109 


-78 


-67 


-70 


-85 


-118 


(5E k + 25m t )/MeV 


+97 


+0.8 


+3.5 


+0.1 


+0.3 












This shows that the prediction for the mass of the IS state is stable in terms of the PS 
mass, and we expect a qualitatively similar conclusion for the top quark cross section. 
(The same observation can also be used to determine the bottom quark MS mass from 
the T resonances accurately |p5[| .) To demonstrate that the PS mass can also be related 
more accurately to the MS mass m t = m^ s (m MS ), we assume fh t = 165 GeV and obtain, 
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adding loop corrections consecutively, 



m t = [165.0 + 7.6 + 1.6 + 0.6 (est.)] GeV (21) 
m tiPS (20 GeV) = [165.0 + 6.7 + 1.2 + 0.3 (est.)] GeV. (22) 



The 2-loop correction follows from |26| and ( |T?D and the 3-loop correction is based on 
an estimate in the so-called 'large- /?oTimit' ||27|| . Hence, the present uncertainty in the 
relation of the PS mass to the MS mass is about 300 MeV. 

4. Discussion. The tt cross section expressed in terms of rrfyps (/•</) with fif — 
20 GeV is shown in the lower panel of Fig. [1]. Since the horizontal scale E is now defined 
as y/s — 2m ti ps(20 GeV), we observe an overall, but trivial shift, related to the fact that 
m t — m t> ps(20GeV) = 1.75 GeV for /i = 30 GeV. The important change is that in the 
PS scheme the location of the peak moves little as we go from LO to NLO to NNLO. 
Furthermore, the scale dependence of the peak location under variations of fi between 
15 and 60 GeV is reduced by a factor of 2. The transition from the pole to the PS 
scheme has little effect on the shape and overall normalization of the cross section as 
expected. In particular a significant uncertainty of about ±20% in the normalization 
remains, larger than at NLO. 

The strong enhancement of the peak for the small scale /i = 15 GeV is a consequence 
of the fact that the perturbative corrections to the residue of the IS pole (see (|T5])) 
become uncontrollable at scales not much smaller than 15 GeV. We find that these large 
corrections are mainly associated with the logarithms that make the coupling run in 
the Coulomb potential. This could be interpreted either as an indication that higher 
order corrections are still important (at such low scales) or that the terms associated 
with bo should be treated exactly, because they are numerically (but not parametrically) 
large. (If the Schrodinger equation with the potential (|II]) is solved exactly by numerical 
methods, the scale dependence of the peak height is indeed smaller [p^fl .) Inspecting 
the logarithms in the result for the cross section, choosing an energy-dependent scale 
/i = 2 (m t (E 2 + I^) 1 / 2 ) 1 / 2 would be most appropriate. Although this choice is our 
preferred one, we refrained from using it for the comparison of the pole and PS scheme, 
since the peak positions are at different energies in the two schemes. 

The uncertainties due to other parameters turn out to be less than the uncertainty 
due to the variation of the scale \i. The dependence on a s {mz) is discussed below. 
Varying fi h by a factor of two around m t changes the cross section by a few percent near 
the peak. The effect of summing logarithms to NLL is of the same order. We have also 
checked the effect of some NNLL logarithms on the calculation and find a variation of 
the order of ±5%. The relatively small effect of renormalization group improvement is a 
consequence of the absence of leading logarithms. We also checked the effect of varying 
[if between 5 and 40 GeV. From a purely pragmatic point of view, values of fif around 
40 GeV lead to the most stable result. However, since /i/ < m t CFOt s {ii) ~ 30 GeV is 
required from a conceptual point of view, we have chosen fif = 20 GeV as our preferred 
setting. (From the point of view of non-perturbative infrared cancellations, it would be 
sufficient to choose /// larger than the strong interaction scale 1 GeV. However, from (|19|) 
we see that the cancellation becomes effective - and the perturbative correction universal 
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Figure 2: Dependence of the NNLO ti cross section on a s (mz) in the PS scheme (solid) and 
the pole scheme (long-short-dashed). The three curves in each scheme refer to a s (mz) = 0.113 
(lower), a s (mz) = 0.118 (middle) and a s (mz) = 0.123 (upper). Recall that E = yfs — 2mt 
in the pole scheme but E = \fs — 2mt t -p$ (20 GeV) in the PS scheme. Other parameters: 
m t = ^ h = 175 GeV, T t = 1.40 GeV, /x = 30 GeV. 

- as soon as the integrand is dominated by \q\ < m t CFOi s . For this reason it is legitimate 
and advantageous to choose \x$ significantly larger than the strong interaction scale.) 

If we take (naively) the change in the peak position under scale variations as a measure 
of the uncertainty of the top mass measurement, we conclude that a determination of 
the PS mass with an error of about 100-150 MeV is possible. Given that the uncertainty 
in relating the PS mass to the MS mass is about 300 MeV, this accuracy seems to be 
sufficient. We emphasize that it is not sufficient to invent an ad hoc mass definition 
in terms of which the peak position is stable empirically. In addition, such a mass 
definition needs to have a well-behaved relation order by order in perturbation theory 
to a mass definition relevant to top quark processes far away from threshold. For a 
realistic assessment of the error in the mass measurement, the theoretical line shape has 
to be folded with initial state radiation, beamstrahlung and beam energy spread effects. 
Since these effects are well understood, the main question that needs to be addressed is 
whether the normalization uncertainty leads to a degradation of the mass measurement 
after these sources of smearing are taken into account. This should be studied in a 
collider design specific setting. 

In the PS scheme the correlation of the top quark mass with a s is also strongly 
reduced, mainly because the perturbative corrections to the IS pole are small in this 
scheme. This is indicated by Fig. ^ which shows the dependence of the line shape on 
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the value of a s (mz)- One can therefore rely less on input from top quark momentum 
distributions, which have been used in the pole scheme to constrain m t and a s (mz) 
simultaneously. Since momentum distributions are more sensitive to non-perturbative 
effects than the inclusive cross section, this is another advantage of the PS scheme. 

In conclusion, we evaluated the next-to-next-to-leading order QCD correction to top 
quark production near threshold in the conventional pole scheme and in the PS scheme 
|1~3|| . We employed factorization in dimensional regularization and summed next-to- 



leading logarithms in the top quark velocity. We find that the cross section expressed 
in the PS scheme allows us to determine the top PS mass more accurately than the top 
pole mass. The PS mass, in turn, can also be related more accurately to the top quark 
MS mass. 
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